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NOMENCLATURE 

C^ Baidwin-Lomax turbulent model parameter 

C kkb Baidwin-Lomax turbulent model parameter 

C, S(dn friction coefficient 

F Baidwin-Lomax outer function 

M Mach number 

u Velocity component in direction of core flow 

u r Friction velocity 

x Location in the direction of the core flow 

y Location in the direction perpendicular to the 

core flow 

v value of y where F is a maximum 

J ma x J 

e eddy viscosity 

ju absolute viscosity 

p density 

x shear stress 

to vorticity 

SUBSCRIPTS 

e core flow value, at the edge of the boundary 

layer 

w values at the wall 

INTRODUCTION 

There is a practical need to model high speed flows 
that exist in jet engine inlets. The boundary layers that 
form in these inlets may be turbulent or laminar and 
either separated or attached. Also, unsteady supersonic 
inlets may be subject to frequent changes in operating 
conditions. Some changes in the operating conditions of 
the inlets may include varying the inlet geometry, bleeds 
and bypasses, and rotating or translating the centerbody. 
In addition, the inlet may be either started or unstarted. 
Therefore, a CFD code, used to model these inlets, may 
have to be run for several different cases. Also, since 


the flow conditions through an unsteady inlet may be 
continually fluctuating, the CFD code which models 
these flows may have to be run over many time steps. 
Therefore, it would be beneficial that the code run 
quickly. Many turbulence models, however, are cum- 
bersome to implement and require a lot of computer 
time to run. since they add to the number of differential 
equations to be solved to model a flow. 

The Baidwin-Lomax turbulence model is a popular 
model. It is an algebraic, eddy viscosity model. The 
Baidwin-Lomax model is used in many CFD codes 
because it is quick and easy to implement. In this pa- 
per. we will discuss implementing the Baidwin-Lomax 
turbulence model for both steady and unsteady com- 
pressible flows. In addition, these flows may be either 
separated or attached. In order to apply this turbulence 
model to flows which may be subjected to these condi- 
tions, certain modifications should be made to the origi- 
nal Baidwin-Lomax model. We will discuss these mod- 
ifications and deter m me whether the Baidwin-Lomax 
model is a viable turbulence model that produces rea- 
sonably accurate results for high speed flows that can be 
found in engine inlets. 

MODIFIED BALDWIN-LOMAX TURBULENCE 
MODEL 

Visbai and Knight [I] suggest several modifications to 
the Baidwin-Lomax turbulence model. Of these, four 
modifications are most significant. These four modifi- 
cations were applied to the "original" Bnklw in-Lomn\ 
turbulence model for tw o dimensional compressible 
steady and unsteady flows that may or may not experi- 
ence separation. These modifications are discussed in 
the following sections. This paper will then evaluate 
these modifications as applied to two cases outlined in a 
later section. 
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Modification 1: and C kW , 

The first modification made to the Baldwin-Lomax 
turbulence model was to the parameters and C kleh . 
The original model proposed by Baldwin and Lomax [2] 
treated these parameters as constants; a value of 1 .6 was 
used for C^, and .3 was used for C^. However, Visbal 
and Knight [l] suggested these parameters are not con- 
stants but are, in fact, dependent on the core flow Mach 
number and may vary by a factor of two over the Mach 
number range 0<M<.<3. In addition, Granville [3] sug- 
gested. if a pressure gradient is present in the flow, 
these parameters should also be modified to include the 
effect of the pressure gradient. He proposed the follow- 
ing relationships for and C^: 

n - 2 - 0.01312 (1) 

klob 3 (0.1724+p) 

V 

q (2 ~4 C kl9b ) ( 2 ) 

( 2 Ckleb~ { 2 ”3 Cfcj •b +( -‘klab ) ) 


where p is the modified Clauser pressure gradient pa- 
rameter defined as: 


0 - 

u T dx 


(3) 


at a streamwise location, x, according to the following 
relation: 


e / -e+<e.-e)e" (jr '*« ,A (4) 


where e 7 is the modified eddy viscosity. \ is the strea- 
mwise location where the pressure gradient becomes 
adverse. is the Baldwin-Lomax eddy viscosity at 
location x a . and A is the relaxation length equal to the 
boundary layer thickness at x,. This relation is valid for 
both the inner and outer regions of the flow. Equation 4 
is used to modify the eddy viscosity when the flow 
experiences an adverse pressure gradient. An additional 
modification was be made to the relaxation length. The 
relaxation length was taken to be a factor times the 
incoming boundary layer thickness at the streamwise 
location just before the pressure gradient becomes ad- 
verse. This multiplying factor is discussed later in this 
paper. 

Modification 3: The Outer Function F 

A third modification that should be made to the Bald- 
win-Lomax model was also suggested by Visbal and 
Knight. They found that the outer function. F. used in 
the Baldwin-Lomax model, was not a suitable length 
scale in regions where the flow was separated or nearly 
separated. This is due to the presence of multiple maxi- 
ma. or peaks, in the F function. The outer function. F. 
is defined as follows: 

F=y | w \D (5) 


The parameter y mnx is used as a length scale, it is the 
location where the Baldwin-Lomax outer function, F, 
has the largest value. This outer function will be de- 
scribed in Modification 3. 

Modification 2: Incorporate Upstream Turbulence 
History Effects 

A second modification made to the Baldwin-Lomax 
model used by Visbal and Knight [I] was to include up- 
stream turbulence history effects, according to a method 
proposed by Shang and Hankey. The eddy viscosity 
was modified in regions with an adverse pressure gradi- 
ent. The Baldwin-Lomax eddy viscosity, s. was relaxed 


where D is the Van Driest damping factor defined as: 

/ ( c\ 

D=1 ' 26 > 


This function. F. will typically have two significant 
maxima in the v-direction. It is necessaiy to choose (he 
correct maximum value of F. F IJ1BX . for each x location. 
Typically, the correct peak is the one furthest awa\ 
from the wall. The value of v ma . chosen, should he ni 
this outer peak. When the boundary layer is not sepa- 
rated or nearly separated, the outer peak will have the 


larger value of F. However, when the boundary layer is 
nearly separated, the inner peak may have the larger 
value of F. However, the outer peak should still be 
chosen for F max and y mv The peak closer to the wall 
will have smaller normal locations, y. Choosing a peak 
that is not the furthest from the wall may yield an 
abrupt and unrealistic reduction in the calculated outer 
eddy viscosity, due to the reduction hi the predicted 

y max ‘ 

There is another problem in determining the maxi- 
mum value of the F function in regions of strong inter- 
action such as separation. In regions of separation and 
downstream of separation, the outermost peak of the F 
function tends to disappear and the innermost peak tends 
to move further away from the wall. This may cause 
large streamwise variations in the normal location of 
F inax , which would in turn yield unrealistic variations in 
the calculated outer eddy viscosity. A way to overcome 
this problem^ to hold the value of the normal location. 
y„, constant for x locations at separation and at all 
streamwise locations downstream of separation. The 
value of y max used should be that of a streamwise loca- 
tion several x stations upstream of separation. F max is 
still calculated for each station downstream of separa- 
tion, however it is calculated at the constant y max (from 
upstream of separation) rather than normal location, y, 
where F is a maximum. 

Modification 4: The Van Driest Damping Factor 

The final modification made to the Baldwin-Lomax 
model was to the Van Driest damping factor. D (equa- 
tion 6). In separated regions, the wall shear stress, t w , 
approaches zero. As Visbal and Knight pointed out, this 
causes an unrealistic reduction in magnitude of the Van 
Driest damping factor, which in turn causes an unrealis- 
tic reduction in the calculated inner eddy viscosity. A 
way to avoid this problem is to use the total shear stress 
at the given normal location, y. to calculate D. This is 
used only for regions where there is separation. The 
total shear stress is determined using the sum of the 
laminar viscosity and the turbulent eddy viscosity. 
However, at a given streamwise station and normal grid 
location, D is itself used to determine the turbulent 
viscosity. Therefore, the known turbulent eddy viscosity 
from the previous streamwise location and same normal 
grid location is used to determine the total shear stress 
to calculate D. 


RESULTS AND DISCUSSION 

The above modifications were made to the Baldwin- 
Lomax turbulence model. These modifications were 
themselves "modified" to better match experimental data 
for high speed flows. Two cases (Table l) were exam- 
ined to obtain numerical results for high speed flows. 
For CASE 1, the Baldwin-Lomax model was incorpo- 
rated in an unsteady boundary layer algorithm that used 
a given inviscid core solution. The Large Perturbation 
Inlet Code (LAPIN), was the code used to determine the 
inviscid core flow solution. The forms of the boundary 
layer equations used in the algorithm were the uncou- 
pled. unsteady, compressible, parabolized Navier-Stokes 
(PNS) equations by Roach, et al [4]. In this case, the 
boundary layer algorithm and LAPIN were interactive. 
The turbulent boundary layer in LAPIN was subjected 
to both favorable and adverse pressure gradients, and on 
occasion experienced separation. Figure 1 shows the 
geometry of the NASA-LeRC 40-60 inlet used for Case 
1. For CASE 2. the modified Baldwin-Lomax model 
was incorporated in an uncoupled, steady, compressible, 
boundary layer algorithm utilizing the PNS equations by 
Roach, et al. In this case the boundary layer was not 
interactive with the core flow. The inviscid core flow 
was prescribed based on the experimental results of 
Lewis, et al [ 5 ]. The turbulent boundary layer in this 
second case experienced no separation but was subjected 
to both adverse and favorable pressure gradients. Figure 
2 shows the core flow Mach number distribution as a 
function of x used for Case 2. In this case, calculations 
were started using a zero pressure gradient profile with 
a Reynolds number, based on momentum thickness, of 
about 4800 at x*11.5 inches downstream of the leading 
edge. 

Table 1 Summary of Cases 
CASE 1 

Unsteady 
Compressible 
Interactive Core Flow 
Core Flow from LAPIN 
Bleeds and Bypasses 

CASE 2 

Stead v 

Compressible 
Non-interactive Core Flow 
Core Flow' Prescribed 
No Bleeds. No Bypasses 
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Effects of Modifying and C Wfb 

In determining the effects of varying C iV and C^ b 
(modification 1) on the Baidwin-Lomax turbulence 
model, modifications 3 and 4 were already incorporated 
into the model. Although modifications 3 and 4 will be 
discussed in greater detail in a later section, it was ne- 
cessary to incorporate them at this point so the empirical 
relations for the parameters and would be com- 
patible with them. 

The original Baidwin-Lomax values of C tT - 1.6 and 
Cueb * -3 predicted C f values that were significantly 
smaller than the experimental C f values for the flow de- 
scribed in Figure 2. (CASE 2). Unfortunately, using 
equations 1-3 in the Baidwin-Lomax turbulence model 
for CASE 2 also predicted C f values that were signifi- 
cantly smaller than the experimental C f values. There- 
fore, another equation needed to be determined to relate 
the parameters C tr and to the core flow Mach 
number and pressure gradient. An equation was ob- 
tained based upon the experimental C f values from 
CASE 2. It was found that the following empirical 
relation: 

1.6 (1+M) (7) 


worked well in flow regions with no pressure gradients 
and flow regions with favorable pressure gradients. 
However, in regions of adverse pressure gradients, using 
the original value of C^* 1.6 compared better to the 
experimental results than using equation 7. Varying the 
constant C Ueb did not significantly effect a change in Q, 
so it was left at its original value of .3. Figure 3 illus- 
trates C f as a function of x for CASE 2 for the follow- 
ing: (1) the experimental results from Lewis et al., (2) 
the turbulence model from Cebeci and Smith [6], (3) the 
Baidwin-Lomax turbulence model with C t7 ,-1.6. and (4) 
the Baidwin-Lomax turbulence model using equation 7 
to determine C^, in regions without an adverse pressure 
gradient. Figures 3 shows that by using the modified 
equation for C^ (equation 7) the model compares well 
to the experimental data and to the model of Cebeci and 
Smith. However the Baidwin-Lomax model, just like 
the Cebeci and Smith model, tends to underpredict the 
skin friction values beyond the adverse region. This is 
consistent with the conclusions made by Visbal and 
Knight. All the turbulence models they tested failed to 
predict the rapid recovery of boundary layers after ad- 
verse regions, such as downstream of the region of 
reattachment after shock waves. 


The above relation for C^ was incorporated into the 
turbulence model for CASE 1, an unsteady, compress- 
ible flow in a NASA-LeRC 40-60 inlet. Figure 4 shows 
a plot of the ratio of the local static pressure to the free 
stream static pressure versus position in the inlet. Using 
Equation 7 to determine C^ as a function of Mach 
number in the favorable pressure gradient region of the 
flow, pushes the shock wave upstream of the shock 
wave position determined from experimental data. Thus 
a further modification to the relationship between C tv 
and Mach number was necessary. This is because Eq- 
uation 7 was empirically developed for a Mach number 
range from 2.5 to 4. The core flow from CASE 1 goes 
from subsonic speeds to about a Mach number of about 
2.5. Equation 5 must be modified to incorporate a Q r 
distribution with Mach number that is compatible for a 
Mach number range from 0 to 2.5. It was determined 
that for lower Mach numbers, below about 2, C^ needed 
to be significantly smaller than the value it had at Mach 
numbers greater than 2.5. The following relationship 
for C tp as a function of Mach number was developed to 
satisfy all these conditions: 

C cp =4 • 6+ ~~ arctan ( 2 . 0 {M-2 .5) ) 

(8) 

Figure 5 shows a plot comparing Equations 7 and 8. 
Equation 8 yields the general trend of C^ as function of 
Mach number for a Mach number range from 0 to 4. 
Equation 7 only captures the general trend for Mach 
numbers greater than 2.5. Figure 4 illustrates that incor- 
porating equation 8 into the Baidwin-Lomax model, 
pushes the shock wave closer to the position determined 
from the experimental data. Equation 8 was also vali- 
dated for CASE 2. Its use in the turbulence model as a 
substitute for equation 7 produced negligible changes in 
the skin friction coefficient. 

It should be noted that if CASE 1 is run for a con- 
stant value of C^ = 1.6. the pressure distribution ob- 
tained is practically identical to the pressure distribution 
obtained when Equation 8 is used to determine C vr . 

This is because the core flow Mach number range is 
from 4 to 2 In this range. C t|l on the most pan averag- 
es a value of about 1.6. thereby giv ing practically the 
same pressure distribution that would be obtained it C ip 
was held constant at the value of 1.6. Equation 8 only 
satisfies the need of transitioning C t1 , to significantly 
higher values, i.e. 4 to 8. for 2.5 < < 4. 
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Effect of Adding Upstream Turbulence History 

Equation 4 was incorporated into the Baldwin-Lomax 
turbulence model for CASE 2. As stated previously 
(modification 2), the relaxation length. A, (of equation 4) 
was multiplied by factor. This factor was varied to 
determine its effect on the skin friction coefficient. The 
factors used were 1, 10, and 25. Figure 6 shows a plot 
of C, versus x for the relaxation length multiplied by the 
factors 1,10, and 25. These results are based upon the 
modified relation for (equation 8). Also included on 
figure 6 are the experimental results of Lewis et aL and 
the results from the numerical scheme of Cebeci and 
Smith. These results indicate that multiplying the relax- 
ation length by 1 yields a skin friction distribution that 
compares very well to the experimental data. However 
using a multiplication factor of 1 yields the same results 
as not incorporating equation 4 in the Baldwin-Lomax 
turbulence model to modify the eddy viscosity in the 
adverse region of the flow. This eddy viscosity modi- 
fication was also incorporated in the region downstream 
the adverse flow region where the pressure gradient 
becomes favorable again for the same cases illustrated 
in figure 6. The turbulence model respectively yielded 
skin friction results that were no different from the 
results illustrated in figure 6. 

The same qualitative results were obtained for CASE 
1 when the relaxation multiplication factors were varied. 
Using multiplication factors of 10 and 25 caused the 
core flow pressure to drop below the experimental val- 
ues in the throat area of the inlet, see figure 7. This 
region is just beyond the adverse flow region which 
modifies the turbulent eddy viscosity according to equa- 
tion 4. Using a multiplication factor of 1 yields a core 
flow pressure distribution in the throat that matches the 
experimental results. This pressure distribution for a 
multiplication factor of 1 is also the same as the pres- 
sure distribution obtained when no modification to the 
eddy viscosity is made (equation 4 is not used). The 
results of CASE 1 and CASE 2 indicate that the use of 
equation 4 is not necessary to obtain results that closely 
correspond to experimental and numerical results. Jn 
fact, use of equation 4 with the relaxation length multi- 
plied by a factor greater than 1 produce results which 
deviate from experimental data. 

Effect of Modifying the Outer Function F 

It was stated in a previous section that the outer func- 
tion, F, had already been modified so as to determine a 
compatible equation for C tT as a function of Mach num- 
ber. This modification insured that the correct maxi- 


mum was chosen at each x- location, in flow regions 
near separation, the outer function, F, may have multiple 
maxima. The correct maximum to choose is furthest 
from the wall. However, in flow regions that are sepa- 
rated and the regions where reattachment occurs, these 
maxima become ambiguous. Calculating F raax and its 
respective y OTILX for x locations at separation and down- 
stream of separation becomes inappropriate. This is 
because F, nax is no longer a valid velocity scale, and y uiax 
is no longer a valid length scale. Visbal and Knight 
determined that downstream of separation the outer 
maxima of F tend to disappear while the inner maxima 
of F, or the maxima of F closer to the wall, move away 
from the surface to a new equilibrium position. Visbal 
and Knight recommend holding y max constant at x loca- 
tions through separation and downstream of separation. 
For CASE I, this constant value of y max was taken to be 
the value of y inax three grid locations, in the streamwise 
direction, upstream of the separation point. Although 
y, ntT was no longer calculated but held at the same value 
for each x station at separation and downstream of 
separation, the outer function F was still calculated In 
this case, the value of the outer function at each x loca- 
tion at separation and downstream of separation was no 
longer the maximum of the outer function. It was sim- 
ply the calculated value of the outer function for the 
given y lUBX and x location. 

It was interesting to see what would happen if for 
CASE I. both F 1MX and y mftX were held constant at sepa- 
ration and downstream of separation. The values of F ttlflX 
and y max would be respectively the values of F um and 
y max three grid locations upstream the separation point. 
Furthermore, it was also interesting to see what would 
happen if neither F loax nor y, mx were held constant at 
separation and downstream of separation, but calculated 
respectively. The results are shown in Figure 8. These 
results indicate that continuing to calculate new values 
for F mM and y mtx at points of separation and downstream 
of separation, positions the shock wave where the exper- 
imental position of the shock wave is. These results 
also correspond very closely to the results of the case 
when only y^ is held constant through separation and 
downstream of separation, while continuing to calculate 
the value of the outer function at this value of y mu> tas 
described above). This case is also illustrated in Figure 
8. Finally, holding F ma , as well as y m . ( at their respec- 
tive values three grid locations in the stieamwise direc- 
tion. upstream separation, pushes the shock wave even 
further downstream of the experimental observation. 
Thus the results from CASE I indicate that there L not 
a significant difference in holding y mas constant tluough 
separation while calculating its respective outer function 
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value. F. and calculating both and through 
separation. 

Effect of modifying the Van Driest Damping Factor 

In regions in and near separation the wall shear stress 
is close to zero, which in turn causes the Van Driest 
damping factor to become small. Consequently there is 
a reduction in the computed eddy viscosity. Visbal and 
Knight recommend using the local shear stress instead 
of the wall shear stress in equation 6 near regions of 
separation and reattachment. 

Figure 9 shows the static pressure distribution versus 
x for CASE I. Only in the region of flow reattachment 
is there a slight difference in using the local shear stress 
as opposed to the wall shear stress in equation 6. For 
the most part the global effects (e.g. core flow pressure, 
shock position) of using the local shear stress instead of 
the wall sheqr stress was not very significant. Using the 
local shear stress in equation 6 was more significant in 
determining skin friction coefficients in region of reat- 
tachment. 

This effect is illustrated in figure 10. Figure 10 is a 
plot of the skin friction coefficient versus x location for 
CASE I. Note that in the region where there is flow 
reattachment, C f is slightly reduced when the wall shear 
stress is used in equation 6. 

CONCLUSIONS 

The Baldwin Lomax turbulence model was modified 
to model unsteady, compressible flows with and without 
separation. Modifications were made to the parameters 
C tp and Cjj^. as well as the calculations of the eddy 
viscosity, the outer function F and the van Driest damp- 
ing factor. The most significant modification made to 
the Baldwin-Lomax model was in making C^ a function 
of both Mach number and pressure gradient. This modi- 
fication was extremely important in matching numerical 
results to experimental data. The other modifications 
did not significantly improve the numerical results as 
did the modification to C . 

The value of C t| , used should be different in regions 
of adverse pressure gradients than regions of favorable 
pressure gradients. In addition, there were two notice- 
able regions of Mach number for using C qr When the 
edge Mach number is below about 2.5. the original 
value of C q = l .6 works fairly well. For edge Mach 
numbers greater than about 2.5. C tp should be signifi- 
cantly larger. Equation 8 seemed to work well for both 


ranges when the pressure gradient was not adverse 
Altering C kleb did not significantly change the results of 
the model, so the original value of C kkh *0.3 was used. 

Modifying the eddy viscosity to account for upstream 
turbulence effects in the regions with adverse pressure 
gradients did not affect the results of the model. Howev- 
er. if the relaxation length scale was increased, the 
predicted results deviated from the experimental results. 

Other modifications to the Baldwin-Lomax turbulence 
model included, holding y mJl constant in regions of 
separation for calculating F^,. and using the local total 
shear stress, instead of the wall shear stress, to calculate 
the van Driest damping factor in regions of separation. 
Both of this modifications improved the accuracy of the 
model. 
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Figure 1 Geometry of the NASA-LeRC 40- 

60 supersonic inlet, used for 
CASE 1. 
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Figure 2 Mach number versus Position down- 

stream, from experiment b> Lewis, 
et al. [5], used for CASE 2. 
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Figure 3 CASE 2: C f versus position for 

C^-1.6, C cp =f(M r , pressure gradi- 
ent) compared with a model of 
Cebeci [6] and experiments of 
Lewis, et al. [5J. 


7 







x/R c 

Figure 4 CASE 1: The ratio of static pres- 

sure to upstream pressure versus 
position for C^'1.6 and calcu- 
lated from equation 7 and from 
equation 8. Experimental results 
from Cubbison, et ah [7]. 
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Figure 6 CASE 2: C f versus position vary- 

ing the relaxation length for up- 
stream turbulence effects. Relax- 
ation length, A, is multiplied by 1, 
10, and 25. 



Figure 5 C^ (vertical axis) as a function of 

Mach number (horizontal axis) 
following equation 7 and equation 
8 . 
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Figure 7 CASE 1: Predicted ratio of local 

static pressure to upstream static 
pressure versus position, varying 
the relaxation length. The relax- 
ation length is multiplied by 1. 10, 
and 25. Experimental observa- 
tions from Cubbison et al. [ 7 ]. 
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Figure 8 CASE 1: Predicted ratio of local 

static pressure to upstream static 
pressure versus position, varying 
the logic of y nmx and F nm calcula- 
tions. Experimental results are 
taken from Cubbison, et al. [7J. 


Figure 10 CASE 1: C r versus position, using 
local and wall shear stress to cal- 
culate van Driest damping factor 
(equation 6). 



Figure 9 CASE 1: Predicted ratio of local 

static pressure to upstream static 
pressure versus position, using 
local and wall shear stress to 
calculate van Driest damping 
factor (equation 6). Experimental 
results from Cubbison et al. [7]. 
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